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We investigate the row sum of the binary pattern generated by the Sierpinski automaton: In- 
terpreted as a time series we calculate the power spectrum of this Sierpinski signal analytically 
and obtain a unique rugged fine structure with underlying power law decay with an exponent of 
approximately 1.15. Despite the simplicity of the model, it can serve as a model for l/f a spectra 
in a certain class of experimental and natural systems like catalytic reactions and mollusc patterns. 
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The phenomenon of 1//" noise is found in a 
widespread variety of systems 0, 0, 0] ■ Usually a noise 
signal is said to be 1// Q (or 1//) if its spectrum follows 
over some decades a power law S(u>) ~ uj~ a with an ex- 
ponent near to one. Despite 1// noise is being known for 
a long time, up to now there is no general explanation 
for universal mechanisms (if they exist). However, on 
a very general level it is believed that complex systems 
are able to generate l/f noise. While real complex 
systems are usually not exactly solvable, extremely sim- 
plified models are under investigation. One prominent 
example is the sandpile-sketching Bak-Tang-Wiesenfeld 
model |5|, a two-dimensional cellular automaton, which 
itself does however not reproduce l/f sufficiently. 

In this work we study an even more simplified model 
introduced in 1984 by Wolfram [|J, which is known to 
be able to exhibit complex behavior: the Sierpinski au- 
tomaton. Looking not at the generated fractal Sierpinski 
gasket itself, but on the row sum, corresponding to the 
total ( in- /activity of the whole system, we have a signal as 
shown in fig Q with increasing mean and increasing spa- 
tial size of the corresponding system; thus every physical 
realization of the system will be finite size-limited. 

Despite on a first glance being a theoretical toy model 
only, Sierpinski patterns have been found in nature. 
Detailed models have explained mollusc patterns b; 
reaction-diffusion models and cellular automata 0, 
Sierpinski patterns also occur in kinkbreedine dynam- 
ics |!j and have been observed in catalysis [l0|. This 
phenomenon occurs generically for suitable parameter 
choices in four standard types of nonlinear spatiotempo- 
ral dynamics including the Bonhoeffer-van der Pol and 
the complex Ginzburg-Landau equation |ll| . 

Consequently, catalytic processes can exhibit similar 
time-signals as the Sierpinski sum signal. A compari- 
son of the reaction rate of a catalytic process with the 
Sierpinski sum signal X(t) = ^Xi{t) has been given in 
[T^ . ll3l |. A single state Xi{t) at a time t can be interpreted 
to indicate the activity of a local catalysis process, i.e. 
reaction (activity) when Xi(t) = and no reaction (in- 
activity) when Xi(t) = 1. The authors |12| observe a 
qualitative similarity between the experimental and the- 
oretical time series. Due to dominating finite size effects 
a 1//" spectra (or long-time correlations) could not be 



identified in the spectrum of the experimental data [Tj| . 
As models of chemical reactions, cellular automata have 
been studied widely |l5j | , explaining spiral waves and pat- 
tern formation in chemical reactions. CO-oxidation on 
Pt(110) and its control by global delayed feedback has 
been studied and compared to models [lfi| , including the 
occurrence of patterns similar to Sierpinski structures in 
the intermittent turbulent phase. Recently, 1/ /" spectra 
have been measured directly in a chemical reaction |l7^ 
by a SQUID setup that allows for much higher resolu- 
tion in time, space and signal-to-noise ratio than the di- 
rect gravimetric measurement of the reaction rate in |l2j ] . 
The power law extends over more than two decades, in- 
dicating the spatiotemporal dynamics of the catalytic re- 
action exhibits avalanches on all sizes and self-organized 
critical behavior . Interestingly, the Sierpinski gasket 
was more recently found in a video feedback system |18| . 
Apart from observation of the Sierpinski pattern itself, 
its geometry has been used widely, e.g., for sandpile dy- 
namics [Tflj l and measurements of magnetoresistance on 
fractal wire networks pcj . 

Definition of the model: The dynamics of the so- 
called Sierpinski automaton is related to the generation 
law of the Pascal triangle. This pattern can be generated 
by the following simple one-dimensional cellular automa- 
ton: We consider a linear array of sites (or spins) Xi(t) 
which can take the values or 1 at discrete time steps t. 
We restrict ourselves to the special initial condition, that 
for t — only one spin is different from all others: 



ar o (0) 



and 



V^o s<(0) = 0. 



(1) 



The dynamics is defined by the following next-neighbor 
interaction: 



Xi (t + 1) = [a*_i(t) + x 4+1 (t)\ mod 2, 



(2) 



i G f— oo, oo], Xi(t) £ {0,1} at dis cre te tir 



time t. In the 

context of catalytic processes [Rj, [12|, 113] a simplified 
chemical interpretation of this rule reads: A catalytic 
process is stopped when too less (i.e. no) or too many 
(i.e. 2) neighbor sites are active. A catalytic process is 
initiated (or continued) when only one neighbor site is 
active. This can origin from a minimal catalysis temper- 
ature combined with a local self-limiting reaction rate. 



2 



The spatiotemporal evolution is obtained from the rows 
of the Pascal triangle by applying modulo 2 elementwisc, 



1 

1 1 
1 2 1 
13 3 1 
1 4 6 4 1 



modulo 2 
> 



1 

1 1 
1 1 
1111 
1 1 



which is also known as the Sierpinski gasket. The Sier- 
pinski gasket is a well-known self-similar structure (with 
point dimension In 3/ In 2 in x for t — > oo) that is ob- 
tained by twofold replication of the first four rows to the 
subsequent four rows, and iteration of this process with 
the whole triangle. 

Instead of considering the fractal pattern itself, we look 
at a scalar observable that can be compared to exper- 
imental time series. Before considering the spectrum, 
we briefly sketch a direct solution |2l| and illustrate the 
analogy to a formal language approach. The row sum (or 
total (in-)activity) over space at time t, defined by 



X(t) :=£>(*). 



(3) 



is referred to as Sierpinski signal. The Sierpinski automa- 
ton rule then generates a time series X(t) (fig-QJ starting 
from t = with 

1, 2, 2,4, 2,4,4,8, 2,4,4,8,4,8,8,16, ... (4) 

Being interested in an analytic expression for if4"jl. we 
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FIG. 1: The self-similar Sierpinski signal X(t) for T = 128. 

first note that X(t) can be generated up to t = 2 N — 1 
from the start sequence uq = (1) by N iterations of the 
sequence replication rule 



u n +i = K, 2 ■ u n ). 



(5) 



Obviously X(t) takes only powers of 2, so we consider 
(lnX(i))/ln2, which starts as 



Strikingly, this appears to be the number of ones in 

0, 01, 10, 11, 100, 101, 110, 111, 1000, 1001, 

1010, 1011, 1100, 1101, 1111, . . . , (7) 

i.e. the binary decomposition of the time variable t start- 
ing with t = 0. Therefore the observable is given by 

X(t) = 2C r ossSum(Binary(i)) (8) 

which is no longer recursively defined. This analytic so- 
lution has already been known in 1852 by Kummer pl| 
in a number-theoretic context. 

A convenient closed expression can be obtained by ex- 
pressing time by a number of spins 

N-l 

t = Y, ° e {o, !}• (9) 

3=0 

and expressing X{{o~i\) from the same configuration as 
X{t) = X(JT a s V\ = 2^=° a \ (10) 

3=0 

By these expressions for i({a"j}) and X({<7j}), we have a 
parametric expression parameterized by a set of N spins 
for all X{t) with times up to i = 2 — 1. 

Spectrum of the Sierpinski signal: The periodogram 
X(u>) of the time signal now is calculated analytically. 
The binary time decomposition allows a Fourier trans- 
formation of X(t) fairly direct from the definition, 

2 N -1 

X(u>) = Y e iut X(t) = Y--J2 e lwt({CT * }) X(*({^})) 



4=0 



N-l 



= II cxp(^(i^+ln2)) 

Co erjv-l j=0 
N-l 

= ]J H eXP (°J ( lUjV + 111 2 )) 

3=0 <?j 
N-l 

= II i 1 + c M^ 3 + In 2)) , (11) 

3=0 

where all sums over Cj are taken over the two possible 
values o~j — and o~j = 1. We now calculate the peri- 
odogram's power spectrum, i. e. |X(cl>)| 2 . The absolute 
value of X(u>) simplifies to a trigonometric product which 
the logarithm converts into a sum, 

N-l 

ln\X(cu)\ 2 = Y ln[5 + 4cos(w2 J ')], (12) 

3=0 

showing a rugged fine structure as shown in fig. [3 A 
rough estimate of the sum in eq. $T2$ is obtained ap- 
proximating the sum by an integral (y :— 

r N-l 

\n\X{uu)\ 2 w / ln[5 + 4cos(w2 J )]dj (13) 



0, 1, 1,2, 1,2,2,3, 1,2,2,3,2,3,3,4, 



(0) 
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ln[5 + 4 cos y] 



dy. (14) 
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As ln(5 + Ax) » ln(5) + |x for |a;| <C 1, we obtain 
, ln5 f 2 ""d« 4 r 2 ""co S (y) 

The integral over the integral cosine is nearly indepen- 
dent of the upper boundary for high values of the bound- 
ary. Thus, we can substitute the upper boundary uj2 n ~ 1 
by some A"-dependent constant, say cm 3> 1. Finally, 
substituting the cosine by one yields immediately a rough 
approximation of eq. I|12") : 



I*MI S 



/ -4/(5 In 2) 



(15) 



Due to the increase of the mean of X(t), spectral estima- 
tion from the periodogram X(uj) (implying a periodical 
extension in time domain) has to be discussed carefully. 

As J2t=o ^(*) — 3 , the average increase is (X(t])L = 
(*(*)}{o...2*-i} - f with /3 = log(3/2)/log(2) » 
0.585. Further, the signal exhibits an increasing variance 
{X{t) 2 ) {0 ... 2N _ 1} = V with 7 = log(5/2)/log(2) « 1.32. 
Consequently, we investigate two variants of X(t): A 
suitable per definitionem mean-free sum signal defined 
by Y(t) — X(t) — (1 + /3)i , and a mean-free signal with 
non-increasing variance 



Z(t) = Y(t)/(Y(t) 



2\l/2 

'{o...t}- 



(16) 



The spectrum of 5^(w) can be directly obtained from eq. 
(|ll|l and the evaluation of the Fourier transform of t@: 



(17) 



The power spectrum of the periodically extended func- 
tion t° decays (for small values of the frequency) like a 
power law with an exponent of approximately -2. Thus, 
the decay is much stronger than X(u>) and the power 
spectra of Y(t) and X(t) deviate only slightly. Hence, 



it follows that \Y(u>) 



\X(u>)\< 



Similarly, 



|X(w)| 2 also estimates |Z(a;)| 2 . We now compare these 
results with numerically applied discrete Fourier trans- 
formation. The power spectrum is fitted (least-squares) 
in fig-Elby a power law with exponent a about 1.11, being 
in good agreement with the analytical result a ~ 1.15 of 
eq. i|15fl • If one measures a power spectrum experimen- 
tally, this may generically be done by observation of reso- 
nances, where the system is coupled with a tuneable oscil- 
lator of given frequency and finite bandwidth. Therefore 
it is quite natural to consider an averaged spectrum. The 
(incommensurable) averaging procedure applied in fig. 01 
smoothes the peaks at ui — 2 k . The peak amplitudes de- 
cay as the average spectrum itself. For commensurable 
averaging (inset of fig. the peaks at ui = 2 k disappear 
completely. 

Note that the spectra \X(cu)\ 2 , and \Y(uj)\ 2 , from FFT 
for T = 2 22 (not shown) also display l/f a behaviour, 
with exponents ax — 1.12, and ay — 1-11 , respectively. 
Moreover, we have used the method of a sliding window 
that normalize the fluctuations |2g of the detrended sig- 

nal Y(t), i.e. Z(t) = Y(t + I) / ' (Y 2 )\ /2 _ l+ht+l} . 
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FIG. 2: Power spectrum of X(t) for T = 1024 time steps 
up to the Nyquist frequency of T/2 (from FFT). The lower 



Envelope is constituted by uj(k) = [- 



,k€ {2,3,...}. 
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FIG. 3: Power spectrum of the time signal Z(t) up to T/8 for 
T = 2 20 time steps (from FFT) and least-square-fit uj' 1 
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FIG. 4: Averaged power spectrum of Z(t) up to T/8 for T = 
2 20 using (incommensurable) l.l fc -bins, i.e. the fc-th interval 
is defined by [[1.1**], [l.l fc+1 ]] where the brackets [*j denote 
rounded integer values. The inset shows the same spectrum, 
averaged using 2 fe -bins, i.e. the fc-th interval is defined by 



[2 fc ,2 



fc+i 



1]. Both correspond to a constant 8uo/u> ratio. 
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For different values of the window width 21, the power 
spectra exhibit power law behaviour with exponents of 
about a » 1.1. Thus 1// Q spectra appear to be a robust 
property of Sicrpinski signals. 

Amplitude distribution: Many systems exhibiting 1/ f 
noise possess a Gaussian amplitude distribution |23. In 
this paragraph we calculate the amplitude distribution 
H N (2 k ) of X(t) analytically where H N (2 k ) denotes the 
frequency occurrence of X — 2 k , k = 0, 1, ... for a sig- 
nal length of T = 2 N . The number of 2 k s in the signal 
up to t — 2 N — 1 is the number of 2 fe s plus the number 
of 2 fc ~ 1 s in the signal sequence up to t — 2 N ~ 1 — 1, i.e. 
H N {2 k ) = H N ^(2 k ) +H N _ 1 (2 k - 1 ). For the boundary 
condition Hn(1) = H\{2) = 1 we obtain a sum of bino- 
mial coefficients H N (2 k ) = £f =1 (I) = * > 1. 



3 

which simplifies for N 3> 1 to 2 N 
distribution: 



N+l\ 

k+l i 

times the Gaussian 



H N {2 k - 1 ) 



(nN/2) 1 / 2 



(k-N/2) z 
N/2 



(18) 



Hence, the amplitude distribution of the occurrence of 
powers of 2 is Gaussian for fixed N. The (averaged) 
amplitude distributions for Y(t) and Z(t) differ from 
i?Ar(2 fc ) but possesses a similar shape as Hx{2 k ). Note 
that the variance distributions for X(t), Y(t), and Z(t), 
are not Gaussian but well defined by eq. ljT%|l . 

As a final point, numerical simulations show that the 
averaged signals are robust against noise, i.e. initial con- 
ditions with more than a single 1. 



In analogous situations in less simply defined systems, 
power laws have also been observed in spatial spectra of 
the scum on fluid surfaces and in the random baker map 
|23l |. and in the temporal spectra in dissipative dynam- 
ics governed by the Lorenz equations |24| being related 
to the Thue-Morse sequence [3 HE l26t l27| . In fact the 
Thue-Morse dynamics 1 — > (1,— 1),— 1 — ► (—1,1) itself 
maps on the string replication rule u n — > (u n , (—1) • u n ) 
for generation of the spatial sequence. Being not equiv- 
alent, but of striking similarity to eq. JSJ for genera- 
tion of the temporal Sierpinski signal series, the Sier- 
pinski dynamics itself does not follow a replication rule. 
While the analytic solution of the Thue-Morse sequence 

is (_1)Cros S Sum(Binary(x)) j^gj ; m ana l ogy to eq . ||, the 

averaged exponents of the resulting spectra are different. 



To conclude, the one-dimensional Sierpinski automa- 
ton generates 1 / f a spectra in the number of active states, 
and can therefore be considered as one of the simplest 
models generating 1// Q spectra. While the Sierpinski au- 
tomaton is rather a caricature, the approach of studying 
the sum signal, or total (in-) activity, and its spectrum, 
can be transferred to more realistic models and com- 
pared directly with experiments. Although exact Sier- 
pinski patterns with long-range correlations remain to 
be experimentally challenging, we conjecture that l/f a 
spectra in a suitable sum signal can be identified in every 
experimental setup exhibiting Sierpinski patterns. 
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